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Abstract 

We analyze the statistics of observables in continuous variable (CV) 
quantum teleportation in the formalism of the characteristic function. 
We derive expressions for average values of output state observables in 
particular cumulants which are additive in terms of the input state and 
the resource of teleportation.. 

Working with a general class of teleportation resources, the Squeezed 
Bell-like states, which may be optimized in a free parameter for better 
teleportation performance [T] we discuss the relation between resources 
optimal for fidelity and for different observable averages. We obtain the 
values of the free parameter of the Squeezed Bell- like states which optimize 
the central momenta and cumulants up to fourth order. For the cumulants 
the distortion between in and out states due to teleportation depends only 
on the resource. We obtain optimal parameters A°|' and Aj'' for the 
second and fourth order cumulants which do not depend on the squeezing 
of the resource. The second order central momenta which is equal to the 
second order cumulants and the photon number average are also optimized 
by the resource with ApJ. 

We show that the optimal fidelity resource which has been found in 
reference [T| to depend also on the characteristics of input tends for high 
squeezing to the resource which optimizes the second order momenta. A 
similar behavior is obtained for the resource which optimizes the photon 
statistics which is treated here using the sum of the squared differences in 
photon probabilities of input and output states as the distortion measure. 
This is interpreted naturally to mean that the distortions associated to 
second order momenta dominates the behavior of the output state for 
large squeezing of the resource. Optimal fidelity resources and optimal 
photon statistics resources are compared and is shown that for mixtures 
of Fock states both resources are equivalent. 

Keywords: Continuous variable quantum teleportation, Photon statistics 



UNIVERSIDAD SIMON BOLIVAR 



1 Introduction 



Continuous variable CV quantum teleportation [2] , depends on the use of a two 
mode quantum entangled resource and in general allows only the reconstruction 
of an imperfect output state. It has been shown that the success probability of 
teleportation can be increased by using entangled non-Gaussian resources which 
can be chosen to improve the efficiency of the protocol [3J 21 [5] • In particular 
for a specific family of two-mode non-Gaussian states, the squeezed Bell-like 
states which include photon-added and photon- subtracted states it has been 
shown that fidelity between input and output states can be improved by carefully 
choosing the parameters of the resource [U [6] . For high squeezing squeezed Bell- 
like states tend to EPR states and teleportation becomes perfect. Other classes 
of non-Gaussian resources, such as two-mode squeezed symmetric superpositions 
of Fock states and of squeezed cat-like states, with good performance in the 
teleportation of single-mode input states were considered in [6j [7] . 

In general, teleportation success is increased by maximizing a functional of 
the density matrices of the output and the input states, such as the fidelity [8], 
over modifications of the resource and the teleportation protocol. A more gen- 
eral approach may be considered since on one side fidelity is not directly observ- 
able except for input states that are eigenstates of an observable quantity and 
in the other any verifying process for the output of teleportation will consist 
in a series of measurements on the output and the input states. Therefore, it 
is worthwhile to study how to optimize the transmission of observable averages 
and/or observable statistics in an output state and to compare such optimal 
transmission situations with optimal fidelity situations. For specific experimen- 
tal situations, the preservation of observable averages and associated variances 
may be in principle of more interest than the overall fidelity. Analysis of the 
optimal teleportation of average values associated with non-classicality for con- 
ditionally produced resources has been done in reference [S], using the transfer 
operator description of CV teleportation [10] . Teleportation of the oscillations 
of photon statistics of squeezed, non-classical states has been done in [11] 

In section ([2]) the characteristic function description of teleportation [12llT|l6] 
is used to obtain averages of observables for the output in terms of contribu- 
tions by the resource and the input. The differences between cumulants and 
raw momenta are stressed. In section ([3j the preservation in teleportation of 
cumulants and momenta is optimized for selected input states over the Squeezed 
Bell-like class of non-Gaussian resources. These optimal preservation resources 
are compared with optimal fidelity resources. In sections Q and ([5]) we use 
the functional T>m of the quadratic deviations of photon probabilities to study 
the distortion of photon statistics through teleportation. We compare T>m with 
Frobenius distance and fidelity. Finally in section ([6]) we present our conclusion. 
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2 Averages of observables for the teleportation 
output 

2.1 The characteristic function description of teleporta- 
tion 

Quantum Teleportation may be described in a very convenient form in terms of 
characteristic functions of the input, the output, and the resource states. For a 
quantum system with complex phase space variable a = x + ip the characteristic 
functions are Fourier transforms of phase space pseudo-distribution functions 
such as Wigner's, P(a) and Q(a) functions. They are generating functions 
for the raw moments associated with their corresponding pseudo-probability 
function [T3]. P(cn) [H] characteristic functions are associated with averages of 
normally ordered products a^ n a m and Wigner's characteristic functions [5J are 
associated with averages of symmetric ordered products. If £ = w + i z is the 
conjugate coordinate to the phase space coordinate a they are related by 

x( e)W = e s i«i 2 / 2 x (o (0) , (1) 

with the symmetric ordering Wigner's characteristic (s = 0) function being 
x(0 = xiO^ an d the P(a) characteristic function (s = 1) for normal ordering 
being x(£) (1) - 

For quantum teleportation in CV, using a symmetric Beam Splitter that 
mixes the A and in states (T=cos(-7r/4) 2 ), with a measurement gain g for the 
homodyne measurement and for the second partner (Bob) correction, any one 
of the characteristic functions for the output state can be written in terms of 
the characteristic function of the input state as, 

Xout (€b ) = xab (gC B ; 6? ) {g£,B ) • (2) 

We use both P(a) and Wigner's characteristic functions in what follows as 
appropriate and take g = 1 for the purposes of this work. 

The output characteristic function in equation @ is a product of the charac- 
teristic function of the input state XinidO an d the characteristic function of the 
resource XAB^g ,£) evaluated at a certain point in the two-mode phase space. 
It is called the Transfer Function of the setup, and has the one-mode phase 
space coordinate £ as argument. For a Gaussian two mode resource state, this 
has been shown to correspond to a one-mode Gaussian characteristic function 
[14], amounting in the output state to a distorting Gaussian noise. 

The distortion of the teleported state in CV depends on the transfer func- 
tion. The transfer function of an state approaching the EPR state is nearly 
constant, being close to the value 1 at each point of the one-mode phase space. 
Thus the limit EPR state is not square- integrable [B]. From equation it is 
seen that with a state approximating an EPR state as a resource the output 
state is nearly identical to the input state. 
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2.2 Averages of operators for the output of teleportation 

For a particular ordering (s) of the creation and annihilation operators a' and 



(a* n a m ) 



(.) 5"+ m xil(0 



(3) 



a + a^)/2 and p = (a — a')/2 i position and momentum 



The averages of x 

operators in CV can be obtained using Wigner's characteristic function, given 
that products x n p m are symmetrically ordered. Thus, with £ = w + i z we 
have, 



(i t n p m ) 



i 



d n+m Xout(0 



(4) 



Now, it is straightforward to calculate the expectation values associated with 
the teleportation output, in terms of the contribution of averages associated 
with the input state, and with the derivatives, evaluated at £ = £* = of the 
transfer function. Though the transfer function may not necessarily be a proper 
characteristic function of a one mode state, for convenience and economy of 
notation their derivatives evaluated at £ = £* = will be denoted as operator 
averages. Given 



(a) n a m ). = 

\ /in 



l 



(x n p m ) 



(i) n + m dC l d£* 



«=e=o 



w—z—0 



(a f n a m ) 



d z n dw m 



d n+m XA B (ge,0 



(x n p m : 



AB 



AB 



i d n+m X AB(ge,0 



n+m 



d z n di 



(5) 
(6) 
(7) 
(8) 



w— z— 



the expectation values for n a m associated with the output state (see equa- 
tions @, ©) are written as, 



i=0,3=0 V / \ ■' / 



AB 



Likewise, for the x n p m averages, we have 



(x n P m ) 



E 

i=0j=0 



m 
j 



-4B 



(9) 



(10) 



The expectation values for the output state are expressed as a sum of prod- 
ucts of " averages" of operators for the input state and the transfer function. 
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For the average of ( n a m ) ou t> we have the terms g n+m ( at n a m )j„ of the 
input state and the average ( at n a m asn + m order terms. The expectation 
value of at n a m for the output state can be expressed as the equivalent average 
for the input state (multiplied by the power of the gain g n+m ), plus terms 
that cause distortion of the input expectation value depending on products of 
lesser-order " averages" of the input state and the transfer function, plus a term 
depending on ( at n a m ) xg. 

The (x n p m ) out expectation value of the output state has a similar form. 
It can be described as the expectation value of the input state (multiplied by 
g n+m ), g n+m (x n p m ) in , plus distorting terms depending on products of lesser- 
order "averages" of the input and transfer function, plus the (x n p m )-jg "aver- 
age" for the transfer function. 

A related, resource-specific analysis for Gaussian resources and photon sub- 
tracted resources of the observable expectation values for the output state of 
teleportation can be found in reference [9]. 

2.3 The transfer function of a two-mode squeezed state 

Wigner's characteristic function of a two-mode resource state prepared from a 
previous symmetric state of modes A and B by two-mode squeezing operation 
with phase of n, can be written as Xab{^'a'^'b) ■> w ^h transformed coordinates 
?a/b(U;Sb) =6i/BCosh(r) -£s /A smh(r). QS]. 

Therefore, the transfer function for such a resource has arguments £a(£* ; — 
£* e~ r and = £ e~ r . It can then be written as 

XAB(?e- r ite- r ) . (11) 

The states prepared by two-mode squeezing have Wigner's characteristic 
functions with the argument |£| 2 . These Wigner's characteristic functions always 
have a multiplying factor of the form e - ^' / 2 , which ensures that they are 
square-integrable, even if they are very non-Gaussian. A transfer function of 
the form outlined above can be written as 

XA B (\ee- 2r ) (12) 

having an envelope of the form e~'^' 2e 2r . The transfer function, being the 
characteristic function of a two mode state evaluated at a particular point in 
two-mode conjugate phase space, satisfies that xab(0',0) = 1 [15] . 

The transfer function is real and symmetric in the conjugate phase space 
defined by £ = w + iz, as exponential and polynomial factors with argument 
|£| 2 e~ 2r are symmetric around £ = 0. While not necessarily Gaussian, the 
transfer function is square-integrable and has a maximum at £ = with value 1; 
given that e - '^ e decreases faster than any polynomial factor in the transfer 
function around £ = 0. This maximum is more salient for small squeezing r. 
For large r, at and near £ = the transfer function is approximately constant. 
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(a) Squeezed Vacuum 



(b) Photon-Subtracted 




(c) Photon-Added (d) Optimal for Coherent 

Input 



Figure 1: The transfer function x(£*>£) with two-mode squeezing r — 1.25 for: 
(a) Squeezed Vacuum; (b) A photon-subtracted resource; (c) A photon-added 
resource; (d) A Squeezed Bell-like resource optimized for the teleportation of a 
coherent input state 

Therefore, the first derivatives of this transfer function, at £ = are 0. 
Second order derivatives are negative or at £ = 0. Thus, for the the transfer 
function of a two-mode squeezed state, 

@)aS = ( & )AB = (*)ab = ®AB = ■ ( 13 ) 

Teleportation using a two-mode squeezed resource does not alter the posi- 
tion and momentum averages (x) in , (p) in of the input state. To illustrate this 
discussion, figure ([1]) displays several sample transfer functions associated with 
well-known teleportation resources. 

The best resources for teleportation are those like the photon-subtracted 
two-mode squeezed state [3] and optimized Squeezed Bell-like state p], which 
have a smoother transfer function around £ = 0. 
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2.4 Output expectation values for relevant observables 
2.4.1 Expectation value of the photon number 

The expectation value of the photon number operator n = a) a can be calculated 
from equation ©, 

(a t a)o«t = g 2 {a) a) in + (a f d)^ 

+ g(a) in (a f >ZB + .9(« f >m ( a >AB ■ ( 14 ) 

The expectation value for h of the input state is modified in the teleportation 
output, by the addition of the "average number of photons" associated with 
the transfer function, and by two cross-terms g (a* ) ( a ) -rg and g (a)i n (a) );j§- 
which are equal to (see equation (1751) ). 



2.4.2 Expectation value of Covariance matrix averages 

The transmission of squeezing by teleportation and the success in teleporta- 
tion of a Gaussian input state may be explored through the comparison of the 
expectation values of the covariance matrix of the output and input states. 

The second-order moments (in x and p, and thus in and a) that make up 
the covariance matrix 

(Ax 2 ) = (x 2 )-(x) 2 (15) 
(Ap 2 ) = (f)-(p) 2 (16) 
Cov(x,p) = (xp) - (x) (p) (17) 

define Gaussian states, up to the average of £ and p or equivalently, up to the 
application of a Glauber displacement operation e a a ~ aa . 
The quotient of the second order central moments 

S = TW)> (18) 

defines the squeezing S of one-mode states. For S ^ 1 states are called squeezed. 
Squeezing of S = e -4 r is produced by applying the transformation 

g2 2 

with a phase 0(£ = r r g R) to a state that initially is not squeezed. 
According to equation (|4]), the averages (a; 2 ) and (x) are given by 

(x)out = g(x) m + (x}^ (19) 
" 2 ^out = g 2 (x 2 ) m + (x 2 )ab + 2g(x) in (x)jg . (20) 



•f- 



A straightforward calculation yields, for ( Ax 



2 



[Ax 2 ) out = g 2 ({x 2 ) m ~ (x) 2 m )+{(x 2 )^ B - (x) 2 -) 

= g 2 (Ax 2 ) m + (Ai 2 )jg , (21) 



G 



where we have defined 



\x)\ 



■ I AB - \ ' AB N ~ ' AB ' 

Equation (f2"Tj) states that the difference between the variances of input and 
output states in CV teleportation does not depend on the teleportation input 
but only on the resource's transfer function. This may be traced to the fact 
that the second order central moment is identical to the second order cumulant. 
The cumulant of order n of a probability distribution that is the convolution 
of two probability distributions, is equal to the sum of the cumulants of order 
n of the two probability distributions [16]. The characteristic function of the 
teleportation output in equation ([2]) is a product of characteristic functions and 
the corresponding Wigner function is the convolution of Wigner functions of the 
input and of the resource [2]. Below, we extend this observation to the third 
and fourth order cumulants. 

The minimization of the distortion of ( Ax 2 ) ou t for any input state can be 
carried out by minimizing the second derivative 



(* 2 )ab = 



d 2 XA B (gC,0 



d; 



w,z—0 



at the origin of the conjugate phase space, since ( x ) 



AB 



An identical argument applies to the variance ( Ap 2 ) out which is given by 
\Af) out = (Af) m + (Ap 2 



i AB 



where, again, we have defined ( Ap 2 



'AB 



'AB 



' AB 



(22) 
For minimum 



distortion in ( Ap 2 } out it is enough to minimize the second derivative 



(f) 



d 2 XA B (gC,0 



AB 



dw 2 



w,z=0 



at the origin of the conjugate phase space. 
Finally, taking into account (|T9]) and 



(xp)out = 9 2 (Xp)i n + [xp /AB 

the average for Cav(x,p) is 



+ 9\x)in (P)ab + 9{P) 



'AB 



(23) 



Cav(x,p)out = g 2 Cov(x,p) in +Cov(x,p) 



AB 



(24) 



where we have defined Cov(x,p)^ = (^p)t§ ~ For the re- 



sources studied is easy to see that (xp\ 



AB 



and in fact Cov(x,p) 



AB 



0. 



There is no distortion of Cov(i, p) ou t produced by a two mode squeezed resource. 



2.4.3 Expectation value of the two-photon correlation function g 2 (0) 

The correlation between the arrival of two successive photons at a photodetector 
for a given time interval is proportional to the media of the product of intensities 
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of the electric field at two different times [T7] . For zero-time delay, this is equal 
to 

5 (0) = 7ataF = 1 + — — (25) 

with Var(n) = ( (a^a) 2 ) — {a) a) 2 being the photon number variance. For clas- 
sical fields, g 2 (0) > f. g 2 (0), is a measure of the non-classicality of the photon 
statistics, fields with non classical photon statistics and <? 2 (0) < 1. ft is worth- 
while to see, then, how g 2 (0) changes in CV teleportation. 

The expectation value {a* 2 a 2 ) can be written using equation ([9]), for a 
normally-ordered output characteristic function: 

(a^ 2 a 2 ) ou \ = g^ 2 a 2 )^ + (a^a 2 n 



AB 



+2. 9 ((at) m (ata 2 )^ + (a) m (at 2 a)W 
+ g ^^ 2 ) m (a 2 ) :rB + ^ 2 ) xs (a 2 ) rn ) 



+2. 9 3 (<at) s (ata 2 ) &> + (a)^(at 2 a)«) . (26) 

The expectation value for (a^a) is given in equation (fT4)l . Note that the first 
order averages (a) jg,(a ') ^jg corresponding to the two- mode squeezed resource 
transfer function (see equation (I13[l ) are equal to zero. In the following section 



v _ AB 1 AB 

two-mode squeezed resource transfer function are also equal to zero. 



we show that the third order averages (a'ci 2 )~ (a' 2 a)^~g corresponding to the 



2.4.4 Expectation value of the central moments of third and fourth 
order 

To study the distortion caused by teleportation in expectation values, for non- 
Gaussian input states, it is necessary to study the expectation values of opera- 
tors which are third or higher order in products of momentum and position. 

The third order central moment of position for the teleportation output is 
given by 

{fJ>£)out = ({X - (£} ou t) 3 ) ou t = (x 3 )out- 3 (x) out (X 2 ) out + 2(x)l u t ■ ( 27 ) 

The averages (x) ou t and (x 2 ) ut have been worked out in equations (|19l) 
and (|20p . From equation (|10[) we obtain, for the third order raw moment 

(x 3 )out = g 3 (x 3 ) ln +3g 2 (x 2 ) m (x) J§ + 3g(x) in (x 2 )^B + (x 3 )^ ■ (28) 

Using the above results, the third order central moment in equation (|2"T|) is 
given by 

(nf)out = (x 3 )^ + 2(x s )^ - Z{ x )~^)~ 
+g 3 ( (x 3 ) tn + 2(x 3 ) ln - 3(x) m (x 2 ) m ) 

= .9 3 <a4 3) >™ + (4 3) )l5 (29) 
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and is an additive quantity of the input characteristic function and the transfer 
function. This is to be expected, since for any probability distribution up to 
third order, all the central moments are identical to cumulants [TB]. An identical 
argument to the one outlined above can be done to obtain, for the third order 
central moment for p 

{nf)out = ff 3 (4 3) )„ + (4 3) )^ • (30) 

The transfer function for a two-mode squeezed resource, as outlined in sec- 
tion (|2.3|) . is of the form x^sd^l 2 e~ 2r ). This function has the third order "raw" 
moment 



dj 2 <97 3 
7 = |e| 2 e- 2 ^ - eHe- 2r = (w 2 +z 2 )e- 2r 



, (31) 

w,z— 



which is easily shown to vanish. Identical reasoning shows that (p 3 )jfg = 0. 
Inspection of equation (|27|) shows that for a transfer function with = 

= 0, the third order central moments (/-4 = (m^)ab = ^- There- 
fore, the third order central moments of input states are not changed by tele- 
portation using a two-mode squeezed resource. 

Furthermore, other third order "raw" moments, such as (x 2 p)^ and (£p 2 )^g 
are also zero. For example 



[X P ' AB " Z3 dz 2 d 1 



(32) 



where 7 is defined in equation (|31j) . Analogous calculations show that (xp 2 )-^ 
vanishes. 

Finally, third order moments of the form (a) " a m ) with n + rn = 3 can be 
written as linear combinations of third order moments of x and p and therefore 
vanish, regardless of operator ordering considerations. Observable quantities 
such as g 2 (Q)out (see equation (|2"5')) '). which are of fourth order in a''' and a, de- 
pend only on the second and fourth order moments associated with the transfer 
function and the input state. 

To further study distortion through teleportation fourth order central mo- 
ments must be analyzed. The fourth-order central moments are not identical to 
cumulants [16] . and hence the fourth order central moment of the teleportation 
output is not equal to the sum of the fourth order central moments of input and 
resource. From 

(A*i ^)out = (( x — ( x ) out) 4 ) out 

= {x A )out - 4 (x) out (x 3 ) out + 6(x 2 ) out (x) 2 out - 3{x) A out , (33) 

(x 4 )out = 9 A {x A ) in + 4g 3 (x 3 ) in (i)jg + 6g 2 (x 2 ) in (x 2 )^ 

+ 4g{x) in {x 3 )^ + (x 4 )^ , (34) 
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it is shown that 



(/4 4) > out = (/4 4) W + (^)ab + % W hn (Ai 2 > 35 (35) 

where (Ax 2 )i n and (Ax 2 )ab are second order central moments. The fourth 
order central moment for p is calculated in an analogous manner; 

(4 4) )o« t = (4 4) ) otlt + (4 4) )j5 + 6. 9 2 (Ap 2 ) m (Ap 2 ) 35 . (36) 

The distortion caused by teleportation on the fourth order central moments 
of the output is a sum of the fourth order central moment of the transfer function 
and a product of the second order central moments of the input and transfer 
function. Thus it is not additive in like order central moments. 

The fourth order cumulant of the output state defined in terms of fourth 
and second order central moments is additive in the corresponding cumulants 
of the input and transfer function [16]. We have, 



and 

As already said 



/v 



4 4) ) ee (4 4) ) - 3 ((Ax 2 )) 2 (37) 
^ = - 3((Ap 2 )f . (38) 



(4 4 )) out = (4 4 ))~ + g\„f) in (39) 
( K f) out = ( K f)~ + g\nf) m . (40) 

For g = 1, the distortion caused by teleportation on the x (p) fourth order 
cumulant of the teleported state is equal to the fourth order cumulant of the 
transfer function. As with the case of lower order cumulants, the distortion 
depends on the resource properties only. Furthermore, the transfer function 
for the two-mode squeezed state described in subsection (12.31) is isotropic in the 



conjugate phase space with coordinate £ = w + 1 z and thus (4^) ab = ( K t^ AB' 
Optimization of teleportation distortion of fourth order cumulants is performed 
through minimization of (4 )ab a l° ne - 

3 Expectation Values Differences and Optimiza- 
tion of Transfer Functions 

We consider now teleportation of chosen input states with the teleportation 
resource given by Squeezed Bell-like states defined by, 

XsuiUl ^)=e- 1 / 2 (l^l 2 + l^l 2 ){A 2 + 2AV1 - A*Re[e i8 &&] 

+ (1-A 2 )(1 - |&| 2 )(1 - |^| 2 )} , (41) 
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with Ca/b = cosn ( r ) — sinh(r) £g , A and fixed squeezing. They depend 

on a free parameter A which can be used to optimize the teleportation output. 
This was done in Q] to optimize fidelity. Here we consider other, more general 
properties of the output and compare the results. 

We consider the following input states: A Fock state with photon number 1, 
a coherent state with a real displacement /3 and a squeezed vacuum state with 
squeezing s. Their characteristic functions are, respectively 

X/iten) = e-*l*»l a (l- I&.I 2 ) (42) 
XcohAZ**) = e-*l*»l a + «M6»]/> (43) 
Xs qc (U) = e-*'«»l a (44) 

£m = fin COsh(s) + Ci„ Smh(s) . 

(45) 

The value of A for optimal fidelity for a Fock state and a coherent state (and 
vacuum state) were computed in [T] for 9 = and g = 1, 

A f" = cos^-arctan^ 3(e 2r_ 1)2 J J (46) 

Kl = cosQarctarxCL + e- 2 '")) . (47) 

The squeezed vacuum optimal fidelity parameter A°p* was also found in Refer- 
ence [1] . 

3.1 Teleportation of Covariance Matrix elements 

Consider the difference between the average ( Ax 2 ) out and ( Air 2 )j n ; 

D A£ 2 = \(Ax 2 ) out - (Ax 2 ) in \ . (48) 

From equation (|2ip. we have 

D AS2 = |(Ai 2 )^| . (49) 

This does not depend on the input. We note that for Squeezed Bell- like 
states the transfer function depends actually on e -2r |£| 2 = e -2r (iu 2 + z 2 ). We 
have (x 2 )^ = (p 2 )jb with 

-)2, 



{S .2y _ d 2 XAB{gC,0 



>AB g z 2 



e- 2r 



2r (6 - 4A 2 - 4Aa/I - A 2 cos(6>)) . 

w,z— 

(50) 

Furthermore, ( x )a~b'(p)a~b as weu as ^ ov (-^p)ab van i sn - We need only to mini- 
mize {x 2 } A ~g which is displayed for a fixed value of r in figure ([2]) This function, 
for 9 = has a minimum at 
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Figure 2: (x )t§ as a function of A and 9 for r = 1.25 

which is independent of the two-mode squeezing r of the resource, although the 
value of (£ 2 )j§ for the Squeezed Bell- like resource goes to as r — > oo. We 
note also that the optimum fidelity resource parameters for the coherent, Fock 
and squeezed vacuum inputs converge towards the same value when r — > oo; 

lim Af* = lim = lim A$ = A°g . (52) 

For large squeezing of the resource the distortion in teleportation of the 
second order momenta dominates the fidelity behavior. 

3.2 Teleportation of Squeezing 

The output state of teleportation has a squeezing equal to 



[Ax 2 ) m + (A£ 2 ) : 
[Ap 2 ) out ~ (A P 2 ) m + (A P 2 /AL: 



<, _ /out _ ivn i / AB . . 



For the two-mode squeezed resources we are considering (Ai 2 )^ = ( Afr ) jg. 
Input states that initially are not squeezed remain unsqueezed and initially 
squeezed input states will become less squeezed after teleportation. 

In figure the quotient S out /Si n between the squeezing of the input and 
output states of a squeezed vacuum state is shown, for a Squeezed Bell-like 
resource with squeezing r = 1.25, and for several values the squeezing parameter 
s of the input state. For the inputs displayed St n — e~ 4s and 1 > S ou t > Si n . 
Note that the quotient is greater for higher values of s. 

The Squeezed Bell- like resource minimizing S ut/Si n is the one that mini- 
mizes the distortion in the output state covariance matrix elements ( Ax 2 ) ou t 
and (Ap 2 ) out , and has A = A°£* ■ 

3.3 Teleportation of Photon Number 

Consider now the photon number average given by equation (| X4|) . For a two- 
mode squeezed resource, "averages" of the transfer function such as (a)^-g = 
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A 

0.75 0.80 0.85 0.90 0.95 1.00 



Figure 3: The quotient of Squeezing S ou t/Si n for a squeezed vacuum state with 
s = (solid line), s = 0.5 (long dashed line), s — 1.0 (dotted line), s = 1.5 (short 
dashed line) and s — 2.0 (dotted-dashed line) as a function of the A parameter of 
the Squeezed Bell-like resource. The two-mode squeezing of the resource has been 
fixed at r — 1.25 



(aT)jg, vanish as discussed in section (|2.3j) . The teleportation distortion of 
photon number average depends only on the resource, regardless of input. The 
optimization of photon number transmission is carried out by minimizing the 
photon number "average" of the transfer function. This is done by minimizing 
the second derivative 



' a* a 



AB 



at the origin of the conjugate phase space. The photon number "average" for 
the transfer function involves a second derivative over £ = w + i z and £* . Thus, 
it is a function of identical behavior to that of equation (|50|) . It is given by 



(fit a) 



AB 



-3 



2 A 2 + 2Ay/ 1 - A 2 cos(0) 



(54) 



and for = 0, it has a minimum at A 



3.4 Teleportation of Fourth Order Central Moments and 
Cumulants 

As discussed above, the fourth order central moments for position and momen- 
tum of the output state are not additive in the central moments of input and 
transfer function. Therefore, the distortion caused by teleportation depends on 
the properties of the input state. Let us measure this distortion as the absolute 
value of the difference between the central momenta, 



D (4) 



(4 



(4)v 



(4 



(4)\ 



(4 



(4)\ 



'AB 



6(Ax 2 ; 



(Ax 2 



AB 



(55) 
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The optimal values of A in this case depend on the two-mode squeezing r 
of the Squeezed Bell-like resource. For the squeezed vacuum input there is a 
dependence on the squeezing s. We have, 



1 



% /(3+e^') 2 (13+2e^(5+ e ^)) 

A ^(4) - J ^ (56) 

(e 2r +3e 2 ») 2 



,/(e 2 ''+3e 2 =) 2 (2e 4 '' + 13e 4 s+10e 2 (''+s>) 

A ^ (4 ) = 1 — 7i ( 57 ) 

3(l+e 2r ) 2 

. + x /(l+e 2 -) 2 (13+30 e 2 - + 18e^) 
A /l,(4) " 7f • 

The squeezed vacuum state is not isotropic in phase-space and has different 
central moments in the powers of x and p. Optimization for the momentum p 
distortion of the fourth order central moment yields 



(3+e 2 (--+°)) 2 

J < 3+e 2 (''+-) ) 2 (l3+2e 2 (''+-) (5+e 2 (-+ = > )) 

A"* = J . (59) 

For s = the equations ([57]) and (|5T)|) reduce to equation ([55]) . The limit for 
r -> oo in equations J5BJ), <J57J) , {53) and is A°£*. 

The transfer function of a non-Gaussian resource such as the squeezed Bell- 
like state has non vanishing fourth order central moments (and cumulants). 
Teleportation using this resource changes the fourth order central moments (and 
cumulants) of any input. Particularly, for a Gaussian input the output will be 
non-Gaussian. On the other hand for g = 1, the teleportation distortion in 
fourth order (x and p) cumulants is determined by ( k ^)xb- For the squeezed 
Bell-like resource with phase 6 = 0, this is given by 



(4 4) )^ = 24e~ ir (-1 + A 2 ) (l - 4Av/l^) (60) 

and is displayed in figure Q. Numerical optimization with respect to the A 
parameter at fixed r yields an optimal value of 

A?£? fa 0.985294 . (61) 

The fourth derivatives with respect to w and z of the transfer function for 
a two-mode squeezed state at w, z = go to zero for large r faster than the 
second derivatives. This explains the convergence of optimal fidelity resources 
when r — > oo to the optimal resources for teleportation of second order moments. 

Cumulants of higher order than fourth go to zero at a faster rate than the 
fourth order cumulant. Therefore, for large r and approximate EPR resources 
the distortion of the teleported state is determined by distortion of second order 
central moments, and resources optimal for second order moments optimize 
fidelity too. 
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Figure 4: The fourth order cumulant («• }xg of the transfer function of a 
Squeezed Bell- like resource with phase = 0, as a function of A and r, 

4 Teleportation of Photon Statistics 

Let us focus now in the differences in the statistics of input and output mea- 
surements for a given observable. A readily measurable quantity in CV is the 
photon number observable, associated with the operator h = a^d, and its eigen- 
states {\n) n = 0, 1, 2, . . . , oo}. 

4.1 Deviations in photon statistics 

Consider the functional 



where the deviations are summed up to an arbitrary finite photon number N, 
which makes T>n numerically computable. For a quantum state described by a 
density operator p and a Wigner characteristic function x(£) the photon number 
probability is given by 



where p n and Xn (£) are respectively, the density operator and the characteristic 
function of the Fock state \n). 

For a quantum state with finite energy (and square-integrable Wigner and 
Wigner characteristic functions) Pn — > when N — > oo. Therefore, the sum 
T>n converges for large N . An upper bound for T>n is obtained when input and 
output are two different Fock states. Then T>n = \/2. 




(62) 




(63) 
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Identical photon statistics corresponding to T>n = does not imply identical 
input and output states, since the phases of the complex amplitudes can be 
different. For example, two coherent states |/3 ) and \/3 e 1 * ) have identical photon 
statistics 

n\ 

though they are different states and can be nearly orthogonal for large |/3| and 
a considerable phase difference <p. A simpler example is constructed with the 
following superpositions of the {|0) , |1)} 

l*+> = ^(|0> + |l» , |*_) = i=(|0> - |l» 

which have Pq = Pi = 1/2 and T> m = 0. 

4.2 Increment and convergence with increasing N 

The increment in T>n associated with an additional ./V + 1-photon term to the 
sum can be approximated by differentiation. Let 

Sn = (P^-P%1{) 2 (64) 

/ N 2 2 \ 2 

\n=0 



= (V 2 N + 5 N y=V N (^ + ^rj ■ (65) 
This can be rewritten as a Taylor series of T>n+i with respect to Sn] 

Vn+1 = Vn + 2^ Sn 82^ + 48^ 4 ~ 384^^ • (66) 

For large enough N, Sjy <C V 2 Nl and quadratic (and higher order) terms in 
equation ([66]) are negligible. Therefore, 

v -»" v " + 4k = v " + 2k «" f - p »" 2 ■ (67) 

We have a positive, decreasing increment for increasing TV. Therefore, T>n 
converges from below. 

4.3 Fidelity and Frobenius Distance 

Fidelity between two quantum states, p m and p out is given by 

F [p m , p*} = Tr (p m p° ut ) = £ E Pnm PZi (68) 



n rn 
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which is equal to the square of the scalar product of the input and output states, 
when both states are pure. It is symmetric to the exchange of density matrices, 
has an upper bound equal to the purity of the purest state involved and a lower 
bound equal to 0. Let p ln = p out , then, 

F[p m ,p m ] = Tr(p m2 ) < 1 . (69) 

Frobenius distance between the same two states is given by 

5 F = (Ti((p out - p m )Hp out - p m ))) 1/2 

= (Tr (p m 2 ) + Tr (p out 2 ) - 2 Tr (p in p out ) ) 1/2 , 

= (Ty(p m2 ) +Ty(p out2 ) - 2F[p in ,p out ]) 1/2 . (70) 

Fidelity can be compared directly to I? at and with the Frobenius distance 
between the two states p m and p out when one of the states (for instance p m ) is 
a mixture of Fock states, with a diagonal density matrix in the Fock basis. 

Let 

p in = J2 p k n) \ k )( k \ • ( 71 ) 

For this state, fidelity is given by 

F[p"\p° ut ] = Y.H P t ] ^\k){k\fn){ m \p out \n) (72) 

n,m k 

= Y, p t ) pV£ = Y, p t ) pt t) ■ (73) 

k k 

Then, 

V N = (P[° ut) f + (Pf l) ) 2 - 2F[/5™,r f ]j ■ (74) 

Using equation ((73)) Frobenius distance can be rewritten; 

5 F = (Tr({p m ) 2 ) +Tr ((p ouf f) - 2 F [p m , p out ] ) 1/2 . (75) 

This is equal to T>n for the special case in which both of the states are 
mixtures of pure Fock states. For inputs which are mixtures of pure Fock 
states, T>n approximates Frobenius distance closely and therefore, V n is a good 
measure of the accuracy of teleportation. 

5 Computation of 

In this section we compute Z?jv of equation (|62[) . up to photon number 24 for 
sample input states and their output states teleported using a Squeezed Bell- 
like resource with varying two-mode squeezing r. The photon statistics for the 
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output state are calculated by numerical integration of the matrix element 



where x«(£) is the characteristic function for a Fock state with photon number 
n. For the input states the exact photon number probabilities are used. 

T>m has been obtained to an accuracy of 7 digits, and found to have a first- 
order differential from 24 to 25 photons that is, for all inputs and for the resource 
squeezing r reported, at least two orders of magnitude smaller than T>n- For 
comparison purposes in the figures, T>n is displayed alongside 1 — F [p m ,p out ] 
(not to scale). 

T>n has been analyzed for Fock states, with photon number (vacuum state) 
and photon number 1 (equation (|4"2"10 ; a mixture of the two Fock states with 
photon numbers and 1 of equal probability P^ n ^ = P± m ^ = 0.5, a coherent 
state (equation (j43|) ) with displacement f3 = 2.12928 and a squeezed vacuum 
(see equation (|44[1 ) with squeezing equal to s = 1.5 

5.1 Vacuum Input State 

The vacuum input state is both a coherent state (with poissonian statistics) of 
displacement 0, and a Fock state with average photon number 0. In figure ([5]), 
T>n and 1 — F are displayed as functions of the Squeezed Bell-like resource's 
A, for four values of the two-mode squeezing r thought to go from a realistic 
resource up to a resource approximating an EPR state. 

It can be seen that the optimal A for D n and 1 — F are noticeably different 
for the smaller values of r. For higher values of r, they become closer. For very 
high r the resource approaches an EPR state, and both optimal As converge 
to A°2*. T> n and the Frobenius distance (see equation ([70]) ) are equal, up to 
the accuracy with which T>n has been calculated. This is in agreement with the 
results in section (|4.3p . for input states that are Fock states. 

5.2 Fock-1 Input State 

The Fock-1 (1 photon number) input state is a highly non-classical state with 
sub-poissonian photon statistics (5^(0) = 0). In figure (j6)), T>n and 1 — F are 
displayed. 

It can be observed that the optimal A for T>m and (1 — F) are different, at 
least for the lower values of r. As r grows, the optimal As converge to A^'. 

As in the case of the vacuum input state, it occurs that T>n and the Frobe- 
nius distance (see equation (JTOj ) are equal, up to the accuracy with which V n 
has been calculated, for the instances of figure (j6|). This is in agreement with 
previous results (see section (|4. 3[) ^) for an input state which is a Fock state. 




(76) 
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(a) r = 0.75 



(b) r = 1.00 




(c) r = 1.25 (d) r = 2.50 



Figure 5: 2?jv (full line, filled circle markers) and (1 — F) (not to vertical scale, 
dashed line, empty diamonds) for a vacuum input state for a Squeezed Bell-like 
resource, as a function of A for: (a) r = 0.75; (b) r = 1.00; (c) r = 1.25; (d) 
r = 2.50 . The central point in each graphic corresponds to the A that maximizes 
teleportation fidelity for the input state , and has been made to coincide in both 
plots 

5.3 Mixture of Fock-1 and Vacuum Input State 

The input state chosen for this example is a mixture of two Fock states, with 
photon numbers and 1, with equal probability P = 0.5. This mixture has the 
characteristic function 

XmUO = 0.5 XcohAC) + 0-5 X /i(0 (77) 

and is a nonclassical state with sub-poissonian statistics (y( 2 )(0) = 0). It is 
also different from all the other input states used in this work in that it is a 
mixed state with a purity (equation (|69p) equal to 0.5. The maximum possible 
fidelity for this input (see section (|4.3|l ) is 0.5. V N and 0.5 — F are displayed in 
figure (JTJ). There is a very noticeable difference, for low r, between the optimal 
A values for T>n and (0.5 — F), which is more pronounced than in the cases 
of the two Fock states previously studied. This difference becomes smaller for 
higher r, as expected, and for very high r with the resource approximating an 
EPR state, there is convergence to A?£) • 

As in the case of the vacuum and Fock-1 input states, T>m and the Frobenius 
distance are equal, up to the accuracy with which T>n has been calculated, for 
the instances of figure (J7]). 
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Figure 6: T>n (full line, filled circle markers) and (1 — F) (not to vertical scale, 
dashed line, empty diamonds) for a Fock-1 input state for a Squeezed Bell- like 
resource, as a function A for: (a) r = 0.75; (b) r — 1.00; (c) r = 1.25; (d) 
r = 2.5 . The central point in each graphic corresponds to the A that maximizes 
teleportation fidelity for the input state and resource used, and has been made 
to coincide in both plots 



5.4 Coherent Input State 

The coherent input state with displacement j3 = 2.12928 has poissonian statis- 
tics, as the vacuum input state, and has the same average photon number value 
(4.534) of a squeezed vacuum with squeezing s = 1.5. In figure (fSJ, Vjy and 
(1 — F) are displayed. 

It can be seen that the optimal values of A for T>m and 1 — F coincide for the 
coherent input state, for the values of r displayed. This is an unexpected result 
for an input state that is neither a Fock state, nor a mixture of Fock states. 
The shared minimum does converge, for higher r, to A"^'- 

The Frobenius distance between the input coherent state and the teleporta- 
tion output is shown in figure (fjJJ). T>n is not equal the Frobenius distance <5p, 
for the coherent input state. However, the optimal A value for the Frobenius 
distance coincides with that for T>n, for all the r values used. 

5.5 Squeezed Vacuum Input State 

The last input state we will consider is the squeezed vacuum input state with 
s = 1.5. It is a highly non-classical state with sub-poissonian photon statistics, 
and has an average photon number equal to that of the coherent state used 
above. T>m and (1 — F) are displayed in figure (fTUj) . 

The results for this input state are very different from those of both the 
vacuum state or the coherent state, for all but the highest values of r. There is 
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(a) r = 0.75 



(b) r = 1.00 




(c) r = 1.25 (d) r = 2.50 



Figure 7: T>n (full line, filled circle markers) and (0.5 — F) (not to vertical scale, 
dashed line, empty diamonds) for a mixture of Fock-1 and vacuum input state 
with equal probability P = 0.5 using a Squeezed Bell-like resource, as a function 
A for: (a) r = 0.75; (b) r = 1.00; (c) r = 1.25; (d) r = 2.5 . The central point in 
each graphic corresponds to the A that maximizes teleportation fidelity for the 
input state and resource used, and has been made to coincide in both plots 



considerable distance between the optimal As for T>n and (1 — F), the greatest 
difference among all the input states used. The expected convergence of A, to 
the value ApJ with increasing r is much slower than for the other input states. 
This is the signature of a highly non-classical input state that is no Fock state, 
nor a mixture of Fock states. 

The Frobenius distance between squeezed vacuum input and teleportation 
output is displayed in figure ([9]). It can be seen that T>n is not equal to the 
Frobenius distance dp. The optimal As for the Frobenius distance, T>n and 
(1 — F) are all different. Only for high r when the resources approach the EPR 
state, do the optimal As converge for the three functionals to A°?. 



21 




(c) r = 1.25 (d) r = 2.50 



Figure 8: T>n f (full line, filled circle markers) and (1 — F) (not to vertical scale, 
dashed line, empty diamonds) for a coherent input state of real displacement 
ft = 2.12928, using a Squeezed Bell-like resource, as a function of the A for: (a) 
r = 0.75; (b) r = 1.00; (c) r = 1.25; (d) r = 2.5 . The central point in each 
graphic corresponds to the A that maximizes teleportation fidelity for the input 
state and resource used, and has been made to coincide for both plots 




(c) r = 1.25 (d) r = 2.50 

Figure 9: Frobenius distance Sp for a coherent input state of real displacement 
/3 = 2.12928, using a Squeezed Bell- like resource, as a function of A for: (a) 
r = 0.75; (b) r = 1.00; (c) r = 1.25; (d) r = 2.5 . The central point in each 
graphic corresponds to the A that maximizes teleportation fidelity for the input 
state and resource used 
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(c) r = 1.25 (d) r = 2.50 



Figure 10: T>m (full line, filled circle markers) and (1 — F) (not to vertical scale, 
dashed line, empty diamonds) for a squeezed vacuum (squeezing s = 1.5) input 
state, using a Squeezed Bell-like resource, as a function of the A parameter of 
the Squeezed Bell-like resource for: (a) r = 0.75; (b) r = 1.00; (c) r — 1.25; (d) 
r = 2.5 . The central point in each graphic corresponds to the A that maximizes 
teleportation fidelity for the input state and resource used, and has been made 
to coincide for both plots 




(c) r = 1.25 (d) r = 2.50 



Figure 11: Frobenius distance 5f for a squeezed vacuum input state of squeezing 
s = 1.5, and the output of teleportation of such an input using a Squeezed Bell- 
like resource, as a function of the A parameter of the Squeezed Bell-like resource. 
The two-mode squeezing r of the resource state is equal to: (a) 0.75; (b) 1.00; 
(c) 1.25; (d) 2.5 . The central point in each graphic corresponds to the A that 
maximizes teleportation fidelity for the input state and resource used 
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6 Discussion 



In this paper we have considered the distortions that teleportation induces on 
second and higher order momenta of gaussian and non gaussian states. Working 
within the characteristic function representation of CV quantum teleportation 
we derive expressions for average values of output state observables in terms of 
the input state and the resource of teleportation. The structure of the output 
characteristic function as a product between the input characteristic function 
and the transfer function and as the Fourier transform of a convolution of the re- 
spective pseudo-distributions leads us to recognize the importance of cumulants 
as the particular functions of raw momenta which arc additive with contributions 
of the input and of the transfer functions. For cumulants distortions produced 
by teleportation are universal in the sense that for any input, they only depend 
on the resource properties. 

We explore teleportation efficiency for cumulants and principal momenta 
using Squeezed-Bell like states as the resource. These states depend on a free 
parameter A and for high squeezing r approximate ideal EPR resources. We 
obtain optimal resources for the different functions and different input states 
and compare with the optimal fidelity resources discussed in Ref. pQ. 

For second order momenta which are also second order cumulants we deter- 
minc from the transfer function of the resource the optimal value A^* = V-Hd^Zl 
independently of the squeezing. For fourth order cumulants we determine nu- 
merically the different optimal value A°^* ss 0.985294. We also determine the 
optimal resources for teleportation of the fourth order central moment for three 
different input states and show that they depend on input state and on the 
squeezing. 

For high squeezing we show that the optimal fidelity resources which also 
depend on the input properties tend to the resource with A"^'. This is explained 
by noting that distortions associated with the higher order momenta decrease 
faster than those associated with the second order moment when r becomes 
large. Even for fourth order momenta which are not cumulants we show that 
the optimal resources for the different input states tend to the one with A°2*. 

We consider the functional T>jq of quadratic deviations and compare photon 
statistics between input and output states. We show that it approximates Frobe- 
nius distance for input states that are mixtures of Fock states. Minimizing T>n 
allows to determine numerically optimal values of A for various inputs. They 
differ from the corresponding maximum fidelity values. As expected, T>n ap- 
proximates Frobenius distance better for Fock states or mixtures of Fock states. 
For every input used and for r arbitrarily large, Frobenius distance, Fidelity 
and T>n optimal values of A converge to A°2*. 

We note that it is straightforward to devise functionals analogous to I? at, 
using other orthonormal, complete bases of eigenstates. These functionals will 
be good approximations to the Frobenius distance for mixtures of states of the 
basis selected. 

Our results support the view that optimal fidelity resources are in first ap- 
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proximation well suited also for an optimal teleportation of more specific char- 
acteristics of the input but that in a realistic situation with r not to large there 
is room for improvement of the resource. 
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